Stochastic Heterostructures in B/N-Doped Carbon Nanotubes 
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Carbon nanotubes are one-dimensional and very narrow. These obvious facts imply that under 
doping with boron and nitrogen, microscopic doping inhomogeneity is much more important than for 
bulk semiconductors. We consider the possibility of exploiting such fluctuations to create interesting 
devices. Using self-consistent tight-binding (SCTB), we study heavily doped highly compensated 
nanotubes, revealing the spontaneous formation of structures resembling chains of random quantum 
dots, or nano-scale diode-like elements in series. We also consider truly isolated impurities, revealing 
simple scaling properties of bound state sizes and energies. 
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The intrinsic transport properties of carbon nanotubes 
are promising for nanoelectronic applications, but use- 
ful intra-tube devices usually require electronic hetero- 
geneity. To date, experimenters have either imposed the 
requisite heterogeneity through electron-beam lithogra- 
phy on a much larger length scale {e.g., gate contacts, 
masked potassium doping (!]]), or exploited preexisting 
inhomogeneities such as on-tube defects Q , chance place- 
ment of external nano-particles, or tube-tube crossings 
jl. One potential route to stable, small-scale hetero- 
geneity of electronic properties is substitutional doping 
by boron and nitrogen. Such doping has already been 
achieved by a variety of techniques, including arc dis- 
charge p|, pyrolysis chemical vapor deposition ||, 
and substitution reactions |7]] . Spatial variation in dopant 
density is much more important in 1-D nanotubes than 
in bulk semiconductors. Fluctuations are large and elec- 
trons cannot bypass anomalous regions. Normally, such 
fluctuations would be considered a nuisance or worse. We 
ask: Can they be exploited? After all, a heavily doped 
compensated nanotube is a chain of p/n junctions (which 
might be called "diodium"). We demonstrate the poten- 
tial of the idea by studying the electronic structure of 
heavily doped (~ 1 atomic %), highly compensated nan- 
otubes. 

This paper takes an atomic-level view of doping and 
junction-like structures. Previously, large-scale inhomo- 
geneous structures, i.e., Schottky barriers and p-n junc- 
tions, were investigated in continuum [jB|-[l0| and self- 
consistent virtual crystal |p"l|| approximations, with an 
emphasis on the meso-scale electrostatic peculiarities of 
one-dimensional systems ||. Such studies provide valu- 
able insights, but cannot access properties arising from 
atomic-scale structure. Ab initio techniques, on the other 
hand, have been applied to a variety of microscopically 
ordered, homogeneous, nanotube materials: BN, BC3, 
and BC 2 N JTJ; BC 40 and NC40 @, C17N Q, and 



abrupt C/BN junctions jlB). There is no previous the- 
oretical study of genuinely randomly (i.e., realistically) 
doped and compensated nanotubes. 

In a (i-dimensional doped semiconductor, the typical 
fluctuation in total dopant charge in a region of linear 
size L is proportional to L d / 2 (assuming Gaussian statis- 
tics), leading to fluctuations in the unscreened Coulomb 
potential proportional to L( d ~ 2 )/ 2 . Only for d = 1 is 
the exponent (d — 2) /2 negative. This implies that the 
dominant fluctuations are then on the shortest length 
scale over which the system is genuinely one dimensional, 
i.e., the circumference. (If the doping is so low that the 
length per dopant [A] exceeds the circumference, then A 
is the dominant scale.) This contrasts with bulk mate- 
rials, for which the most important fluctuations can be 
on very long length scales and are controlled by screen- 
ing, especially for light doping and high compensation 
]l6[ . Screening is ineffective at short distances in nan- 
otubes; taking screening into account merely strengthens 
the conclusion. We can estimate the potential fluctua- 
tion (AV) over an axial distance L = 2~kR by computing 
the average potential produced by a Gaussian ("%/iV") 
charge fluctuation. Denoting the proportion of atoms 
which are dopants by c, the dielectric constant by e [fl7| , 
and using the graphene lattice constant 2. 45 A, we find 
AV ~ 9 eV c 1 / 2 /e. This is consistent with the numerical 
results presented below. Note that AV is independent of 
R for fixed c. Fluctuations in the azimuthal locations of 
dopants are not important: all azimuthal Fourier com- 
ponents of the doping have the same variance, but the 
potential arising from nonuniform components (m 7^ 0, 
see Eq. |l|) falls off rapidly and cannot be resolved by 
intra-band processes. (Numerical results confirm that 
they have little effect.) 

Ab initio methods are prohibitively expensive for large 
systems, and ordinary tight binding is unable to treat 
charge transfer over distances greater than about a bond 
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FIG. 1. Test of SCTB against LDA on a regular planar 
BC10N structure. Top: the unit cell. Bottom: Electronic 
energy bands near the Fermi level (dashed line) calculated by 
simple TB, SCTB, and LDA-DFT. 



length, so that we have chosen to use a self-consistent 
tight-binding (SCTB) technique. The nonorthogonal fl§| l 
tight-binding basis contains four orbitals per atom with 
Hamiltonian and overlap parametrization determined di- 
rectly from LDA @. This parametrization was shown 
to give results in good agreement with LDA over a wide 
range of structures with short-range charge transfer. We 
impose self-consistency of long-wavelength components 
of the Hartree Coulomb potential via an iterative pro- 
cedure |^| . Figure [j] compares results of simple TB, 
self-consistent TB, and LDA for a strongly charge mod- 
ulated planar structure. SCTB is remarkably effective 
at correcting the deficiencies of simple TB, reproducing 
the main band-narrowing and shift of Fermi energy level 
crossing exhibited by the LDA calculation, and is there- 
fore expected to describe the physics near E F well. The 
data presented here are for mechanically unrelaxed sys- 
tems. Tests show that relaxation introduces only minor 
corrections to the electronic structure. 

A length of nanotube heavily doped with both boron 
and nitrogen might have a substitution pattern like that 
shown in the second panel of Figure g. This shows a sec- 
tion of (8,0) nanotube (cut open and rolled flat with the 
axial direction horizontal) containing a number of boron 
and nitrogen atoms substituted on randomly chosen sites. 
The section is treated as the unit cell of a superlattice. A 
more lightly doped sample is shown in figure |^. Fluctu- 
ations in the doping act as collective impurities, creating 
localized impurity states in the energy gap of the un- 
doped nanotube. This localization is illustrated by the 
plots of the axial probability densities (azimuthal aver- 
age of |V>| 2 ) of deep collective donor and acceptor states. 
These states arise from a local imbalance in the dopant 
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FIG. 2. SCTB results for sample A, an (8,0) nanotube ran- 
domly doped with 10 each nitrogen and boron atoms per su- 
percell of length 7.7 nm (overall doping level 3.5%). The 
horizontal direction in all figures corresponds to the axial di- 
rection, (a) Azimuthally averaged self-consistent potential. 

(b) The doping pattern in the unit cell cut open and laid fiat. 

(c) Axial probability density profiles of deep donor and ac- 
ceptor states (marked by arrows in the DOS plots) . (d) Total 
densities of states for the corresponding slices indicated by 
dashed lines in panel (b). Dashed lines show the gap edges of 
an undoped (8,0) nanotube (centered at zero). 



charge of only one or two unit charges. Regions with ex- 
cess nitrogen atoms suppress the local density of states 
near the bottom of the gap {vice versa for boron-rich 
regions). This remnant of band bending would matter 
little in a bulk material, but in a one-dimensional struc- 
ture can effectively block charge transport if the Fermi 
level lies in such a local pseudogap. 

These randomly doped nanotubes resemble disordered 
chains of tiny quantum wells or dots, and could show 
a variety of behaviors characteristic of such structures, 
e.g., negative differential resistance. Level spacing in the 
wells is A ~ lO^-Egap. Band widths for the well states, 
r ~ 1 meV for the samples studied, measure interwell 
tunneling matrix elements. The charging energy scale, 
Ec = e 2 /(e2nR) — (0.64/e)£ , gap , scales with the gap. 
This is a significant energy, so that one must remember 
the local DOS plots in the Figures are specific for neutral- 
ity in the depicted unit cells. A length of compensated 
doped semiconducting tube, particularly under imper- 
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FIG. 3. Sample B. SCTB results for an (8,0) nanotube with 
a better segregated and lower density arrangement than sam- 
ple A, containing 4 each nitrogen and boron atoms per 7.7 
nm unit cell. Presentation is as in Fig. El 



feet dopant mixing, forms a dense sequence of alternately 
directed nanoscale diode-like elements. Reverse-biased 
structures should dominate transport under a voltage; 
field doping the tube would wash out weaker junctions, 
isolating regions of most pronounced nonlinear behavior. 

Leonard and Tersoff j8| provided a formula showing 
that the depletion layer of a nanotube p/n junction is 
very long at low doping, but shrinks exponentially with 
doping. Unfortunately, the formula is not applicable to 
the high doping (~ 1%) regime relevant to devices with 
nanometer-scale structure. The depletion width in sam- 
ple B is close to being saturated and is very narrow (~ 1 
nm). 

Remarkably, many features of this picture also hold 
for doped metallic nanotubes, as shown in Figure || The 
Fermi level of a pristine metallic tube lies at the center of 
a subband which contributes a small and constant density 
of states. Upon doping, this subband forms a relatively 
inert background over which appears a response very sim- 
ilar to that of a semiconducting tube. One observes band 
bending and highly concentrated (though not technically 
localized) states in the would-be gap. The roughening 
of the local DOS near the Fermi level seen in the figure 
is partly an artifact of periodicity and k-point sampling, 
but it also suggests that backscattering may be stronger 



than in an extrapolation from the limit of vanishing dop- 
ing@. 

We lastly consider states bound to isolated impuri- 
ties, using both SCTB and an effective mass approxi- 
mation. These studies reveal new diameter-dependent 
scaling relations j22]. A surface charge density p m — 
(Qm/27r)<5(z)e™^ on a cylinder of radius R (cf> is the az- 
imuthal angle, z, axial position) produces a potential at 
(z, (j) = 0) which falls off asymptotically (z — > oo) as 
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This familiar result of multipole expansion shows that 
only the m = component of charge has an influence of 
any appreciable range. We retain only the m = com- 
ponent of the impurity Coulomb field, partly due to this 
estimate, but mostly because m ^= contributes only 
if there is band mixing. The m — component of the 
potential can be written explicitly as an integral, but it 
differs little from the potential of an identical charge lo- 
cated on the tube axis, so we use the latter for simplicity, 
taking the charge to be Z\e\. Ignoring the trigonal distor- 
tion of the nanotube electronic structure, the band-slicing 
picture gives an effective mass m* — E gap / (2h 2 v 2 ,), with 
which we obtain the one-dimensional wave equation 
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where £ = z/R is the axial coordinate scaled by the tube 
radius. The bandgap of a large-gap nanotube (one whose 
wrapping indices do not differ by a multiple of three) 
satisfies E gap R — 2hv F /3. Therefore, if Z/e is held fixed, 
both E/E mp , the ratio of bound state energy to the gap, 
and ip(z/ R) are independent of tube radius. Surprisingly, 
the bound state size and energy scale with the diameter 
of the tube. 

The energies and axial spreads of the bound states are 
determined numerically from Eq. @. Empirically, both 
scale extremely well with Z/e between 1 and 10. We find 



~\1.25 



E/E gap » -2M(Z/e) 
((z/R) 2 ) « 1.2(Z/e)~^ 
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These results are consistent with self-consistent tight- 
binding. Our SCTB calculation for an (8,0) nanotube 
with one nitrogen or boron impurity per 7.7 nm unit cell 
verifies the existence of a bound state fairly well localized 
within that region ((Az) ~ 3.2R). Since the calculation 
assumes a periodicity, the average potential in the unit 
cell does not correspond to that of an isolated impurity, 
so we must subtract the relevant offset. This procedure 
is somewhat uncertain since e(k) varies from 7 to 4.5 as 
kR falls from 1.0 to 0.3. Assuming k — 27r/(cell length) 
gives a binding energy of about 0.23E g:ip . This value is 
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FIG. 4. Sample C. SCTB results for a doped (6,6) nan- 
otube (6.9 nm unit cell) with 1.2% each boron and nitrogen. 
The dashed lines, centered at zero, indicate the edges of the 
m = ±1 subbands of the undoped tube, at which the DOS 
jumps suddenly from a low constant value. 

duplicated by Eq. (||) with 1=1, which is reasonable due 
to the importance of the short-range part of the poten- 
tial. With e — 7, the model predicts a spread (Az) ~ 3R, 
very close to the SCTB result. 

In summary, the pointlike character of dopant charges 
in B/N substitutionally doped nanotubes is important 
at all doping levels. At high doping, fluctuations sponta- 
neously produce a high density of nonlinear nanometer- 
scale quantum dot and/or junction structures, which are 
relatively insensitive to environmental perturbations such 
as packing into nanotube bundles (note that these junc- 
tions do not simply provide scaled-down traditional func- 
tionality) . These results are also relevant to chemisorbed 
dopants. The one-dimensional nature, which makes dop- 
ing variation especially relevant to electronic character- 
istics, also facilitates mapping and manipulation, e.g., 
by scanning tunneling microscopy. Using fluctuations in 
this way can produce nonlinear device characteristics on 
scales far smaller than can be fabricated intentionally and 
avoids the need for fine doping control. 
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